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ABSTRACT 

Photometry has been obtained for magnitude limited samples of galaxies in the two 
rich clusters Abell 665 (37 galaxies) and Abell 2218 (61 galaxies). Both clusters have 
a redshift of 0.18. The limiting magnitude of the samples is 19™ in the I-band. Spec- 
troscopy has been obtained for seven galaxies in A665 and nine galaxies in A2218, 
all of which also have available photometry. Spectroscopy has been obtained for two 
additional galaxies in A2218, one of which is a background galaxy. 

Effective radii and mean surface brightnesses </>c were derived from the pho- 
tometry. The typical uncertainties are ±0.078 in logrc and ±0.12 in log</>c. The 
combination log To + 0.82 log </>c that enters the Fundamental Plane (FP) has an un- 
certainty of only ±0.018. The spectroscopy was used for measurements of the velocity 
dispersions, a. The typical uncertainty is ±0.023 on logcr. 

The data are used to establish the FP, log To = a logcr + /31og</>e ± 7, for the 
clusters. The FP for these two clusters adds important knowledge about the properties 
of E and SO galaxies in the relatively unexplored redshift interval 0.05 to 0.3. We have 
compared the FP for A665 and A2218 with the FP for the three clusters CL0024+16, 
CL1358+62 and MS2053-04 with redshifts between 0.33 and 0.58, and with the FP 
for the Coma cluster. The scatter around the FP is similar for all six clusters. We 
find that the FP for the intermediate redshift clusters has a smaller coefficient a than 
found for the Coma cluster and other nearby clusters. This may either be caused by 
selection effects for the intermediate redshift clusters or by differences in the evolution 
of low luminosity galaxies and high luminosity galaxies. 

The mass-to-light (M/L) ratios, as measured from the FP, change with redshift. 
At z=0.18 the M/L ratio for photometry in Gunn r in the rest frame of the clusters 
is (11 ± 7)% smaller than for the Coma cluster. Using the data for A665 and A2218 
together with the previously published data for CL0024+16, CL1358+62 and MS2053- 
04, we find that the M/L ratios for photometry calibrated to Gunn r change with 
redshift as A log M/L^ = (-0.26 ± 0.06) Az for = 0.5. This change in the M/L ratio 
is equivalent to the absolute magnitudes changing as AMj-.^. = (— 0.65±0.15)Az. These 
new results are consistent with the previously published analysis for CL0024+16, 
CL1358+62 and MS2053-04. 

For (7o = 0.5 the results are consistent with passive evolution of a stellar population 
which formed at a redshift larger than five. For Qo = 0.15 the formation redshift must 
be larger than 1.7. Our new data for A665 and A2218 confirm the gradual and slow 
evolution of the bright E and SO galaxies. However, possible star formation in E and SO 
galaxies during the last 3-4 Gyr of the history of the Universe, as well as the selection 
effects complicate the interpretation of the data. 

Key words: galaxies: clusters: individual: Abell 665, Abell 2218 ~ galaxies: elliptical 
and lenticular, cD - galaxies: evolution - galaxies: stellar content. 
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1 INTRODUCTION 

Observational studies show that the formation and evolu- 
tion of galaxies is a complex process, which may involve in- 
teractions, star bursts and infall (e.g.. Dressier et al. 1994ab; 
Lilly et al. 1996; Moore et al. 1996). It has also been found 
that some nearby E and SO galaxies may have experienced 
star formation in the last 3-4 Gyr. Caldwell et al. (1993) 
found some E and SO galaxies in the Coma cluster to have 
post-starburst spectra, and Faber et al. (1995) suggest that 
nearby field E galaxies have a substantial variation in the 
mean age of their stellar populations. 

In order to study the evolution of galaxies we need to 
study both the morphological evolution as well as the evolu- 
tion of the luminosities and the mass-to- light (M/L) ratios of 
the galaxies. High-resolution imaging with the Hubble Space 
Telescope (HST) and from the ground, combined with spec- 
troscopy from the ground make it possible to carry out this 
kind of studies. Studies of the morphological evolution show 
that the fraction of spiral galaxies and irregular galaxies in 
clusters was higher at larger redshifts (e.g.. Dressier et al. 
1994ab; Oemler et al. 1997; Dressier et al. 1997; Couch et 
al. 1998). A detailed investigation of how the luminosities 
of disk galaxies change with redshift is possible using the 
TuUy-Fisher (1977) relation. Vogt et al. (1996) established 
the TuUy-Fisher relation for a sample of field galaxies with 
redshifts between 0.1 and 1 and found a luminosity evolution 
in the B-band of « 0^6 between z — 1 and the present. 

The Fundamental Plane (FP) (Dressier et al. 1987; 
Djorgovski & Davis 1987) for elliptical galaxies makes it 
possible to study how the luminosities and the M/L ra- 
tios of these galaxies change with redshift. Also, SO galaxies 
in nearby clusters follow the FP (e.g., J0rgensen, Franx & 
Kjffirgaard 1996, hereafter JFK96). The FP relates the ef- 
fective radius, r^, the mean surface brightness within this 
radius, </>e, and the (central) velocity dispersion, a, in 
a tight relation, which is linear in log-space. For 226 E 
and SO galaxies in nearby clusters JFK96 found logrc — 
1.24 logo- — 0.82 log <7>o -I- est, for photometry in Gunn r. 
This relation may be interpreted as M/L oc M°-^*rc"° °^ 
(cf., Faber et al. 1987). The change of the M/L ratio with 
the mass is partly caused by changes in the stellar popula- 
tions. Some of the apparent change may also be explained 
by non-homology of the galaxies (e.g., Hjorth & Madsen 
1995). Throughout this paper we assume that the FP im- 
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plies a relation between the M/L ratios and the masses of 
the galaxies. 

The scatter of the FP is very low, equivalent to a scatter 
of 23% in the M/L ratio (e.g., JFK96). Thus, the FP offers 
the possibility of detecting even small differences in the M/L 
ratios by observing a handful of galaxies in a distant cluster. 

The FP has been used to study the evolution of the 
M/L ratios of cluster galaxies as a function of redshift up 
to a redshift of 0.83 (van Dokkum & Franx 1996; Kelson 
et al. 1997; Bender et al. 1998; van Dokkum et al. 1998; 
Pahre, Djorgovski & de Carvalho 1999). All of these studies 
find the M/L ratios of the galaxies to increase slowly with 
decreasing redshift. Only the study by Pahre et al. contains 
clusters with redshifts between 0.05 and 0.3. In this paper 
we establish the FP for the two rich clusters Abell 665 and 
Abell 2218. Both clusters have a redshift of 0.18, adding 
data to a redshift interval not well-covered by the previous 
studies. 

We use the data for A665 and A2218, data for the three 
clusters studied by van Dokkum & Franx (1996) and Kelson 
et al. (1997), and data for the Coma and the Hydral clusters 
to study how the FP and the M/L ratios of the E and SO 
galaxies change with redshift. The full sample covers about 
half the age of the Universe. 

The available photometry and spectroscopy, as well as 
the derivation of the photometric and spectroscopic parame- 
ters is described in Sections 2 and 3, respectively. The details 
of the data reduction are provided in Appendix A and B. 
Appendix A also contains photometric parameters (effective 
radii, mean surface brightnesses and colors) of magnitude 
limited samples of galaxies in A665 and A2218. Section 4 
briefiy describes the data for the other clusters. The FP is 
analyzed in Section 5. In this section we also discuss the 
evolution of the M/L ratios of the galaxies. The conclusions 
are summarized in Section 6. 



2 THE PHOTOMETRY 

The central parts of the two rich clusters Abell 665 and 
Abell 2218 were observed with the Nordic Optical Telescope 
(NOT), La Palma, in March 1994. Observations were done in 
the I-band and the V-band. The sizes of the observed fields 
are 3.'4 x 2!7 for A665 and 2.'8 x 3!5 for A2218. Table § 
summarizes the instrumentation for the observations. The 
I-band images were taken as individual 600 sec exposures, 
with a few shorter exposures. The total exposure time, the 
seeing and the sky brightness of each of the final co-added 
images are given in Table ^. The V-band images were taken 
on another night than the I-band images, and have rather 
poor seeing. 

A665 was observed with the HST with the Wide-Field- 
Planetary-Camera 2 (WFPC2) in the fihers F606W and 
F814W. Two fields were observed giving a total coverage 
of approximately 9.4 square arc-minutes on the WF chips. 
For A2218 we use HST/WFPC2 archival data obtained in 
the F702W fiher. One field in Abell 2218 was observed with 
the HST, covering approximately 4.7 square arc-minutes on 
the WF chips. The images were taken as individual expo- 
sures with exposure times between 1600 sec and 2300 sec. 
The total exposure times for the co-added images are given 
in Table |[ 
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Table 1. Instrumentation 



Photometry Spectroscopy 



Dates 


March 4-9, 1994 


Sept.-Dec, 1994 


Jan. 18-19, 1991 


June 6-8, 1994 


Telescope 


NOT, LaPalma 


HST 


MMT 


KPNO 4m 


Camera/Spec. 


Astromed 


WFPC2 


Red Channel Spec. 


R&C Spec. 


CCD 


EEV 1152x770 


Loral 


TI5 


T2KB 


r.o.n. 


4.5 e- 


5.0 e- 


9.0 e- 


7.4 e- 


gain 


0.9 e-/ADU 


7.5 e-/ADU 


2.0 e-/ADU 


2.6 e-/ADU 


Pixel scale 


0.1635" /pixel 


0.0996" /pixel" 


0.6"/pixel 


0.69" /pixel 


Slit width 






1.8" 


1.9" 


Spectral res., cr 






1.7lA 


1.38A 


Wavelength range 






5780-6420A 


4255-6970A'' 


Cluster (s) 


A665, A2218 


A665, A2218 


A665 


A2218 



Notes. ° All galaxies studied in this paper are imaged on the WF chips of WFPC2. ^ The unvignetted 
range is 4550-6500A. 



Table 2. NOT 

Image 


images 
Filter 


^cxp 

[sec] 


seeing 
["] 


"Isky 

[™/arcsec^] 


A665 field 1 


I 


9600 


0.66 


20.0 


A665 field 1 


V 


1200 


1.86 


21.8 


A665 field 2 


I 


5874 


0.60 


19.6 


A665 field 2 


V 


1200 


2.03 


21.8 


A2218 field 1 


I 


5400 


0.70 


19.3 


A2218 field 1 


V 


1200 


3.33 


21.2 


A2218 field 2 


I 


5400 


0.82 


19.4 


A2218 field 2 


V 


1200 


4.23 


21.2 



Notes. The seeing is measured as the full-width-half-maximum 
(FWHM) of a Gaussian fit to the core of images of stars. A665 
field 1 is centered 0.'6 east and 0.'3 south of the BCG; field 2 is 
oSset 1.'7 to the west of field 1. A2218 field 1 is centered 0.'3 to 
the south of the BCG; field 2 is offset 1.'7 to the north of field 1. 



Table 3. HST 


images 




Image 


Filter 


^cxp [sec] 


A665 field 1 


F606W 


4400 


A665 field 1 


F814W 


4400 


A665 field 2 


F606W 


5100 


A665 field 2 


F814W 


4800 


A2218 


F702W 


6500 



Notes. A665 field 1 is centered 0.'3 south-east of the BCG. A665 
field 2 is centered 1.'7 west of the BCG. The observation of A2218 
is centered 0.'9 east of the BCG. 

The basic reductions of the images are described in Ap- 
pendix A, which also covers the calibration of the photom- 
etry to the standard passbands 7c and V. In the following 
we describe the determination of the global photometric pa- 
rameters (effective radii and mean surface brightnesses) and 
the calibration of the photometry to Gunn r and Johnson 
V in the rest frames of the clusters. 



2.1 Global parameters 

Determination of the effective radius, r^, and the mean sur- 
face brightness inside this radius, <ij.>^, was done following 
the technique described by van Dokkum & Franx (1996). 
This technique uses the full 2-dimensional information in 
the image. Each galaxy is fit with a 2-dimensional model. 



taking the point-spread- function (PSF) into account. The 
output from the fit is the center of the galaxy {x,y), the 
effective radius, r^, the surface brightness at the effective 
radius, Je, the ellipticity, e, the position angle of the major 
axis and the level of the sky background. All galaxies pre- 
sented here were fitted with r^/** profiles. The mean surface 
brightness <IJ,>^ was derived from the local surface bright- 
ness at Tc under the assumption that the galeixies have r^^* 
profiles. 

The PSF images for the NOT I-band images were con- 
structed from 2-3 well-exposed stars in each image. The PSF 
images have a size of IC/.'S x lO'.'S. It was tested that the PSFs 
did not vary over the field. The PSF images for the HST ob- 
servations were constructed using Tiny Tim (Krist & Hook 
1997) and have a size of 3" x 3". 

All galaxies brighter than an apparent total magnitude 
in the I-band of mx ~ 20™ were fitted. The brightest cluster 
galaxies (BCG) in both clusters were fitted within a radius 
of l-2ro. The fainter galaxies were fitted within the radius 
where the local surface brightness reached ~ 26™/arcsec^ in 
the I-band. This radius is typically 5-7rc. 

Before fitting the fainter galaxies the BCGs were mod- 
eled with the program Galphot (Franx, lUingworth & 
Heckman 1989a; j0rgensen, Franx & Kjaergaard 1992). 
This program fits a model of ellipses to the isophotes of 
the galaxy. The model was then subtracted. The use of 
Galphot ensures that gradients in the background due to 
signal from the BCG are properly removed before fitting 
the fainter galaxies. Further, the isophotes for the BCG in 
A2218 are not co-centric, as also noticed by Kneib et al. 
(1995). Galphot handles this properly while the fitting 
technique described above uses a fixed center and there- 
fore leaves large residuals, which in turn affect the results 
for the fainter galaxies. The BCGs were also modeled with 
the 2-dimensional fitting technique. The effective radii de- 
rived from this agree within 0.02 in logTe with the results 
from fits to the growth curves produced by Galphot. Fit- 
ting the growrth curves produced by Galphot does not give 
the possibility of 2-dimensional modeling the effect of the 
PSF. In the following, we therefore use results from the 2- 
dimensional fitting technique. 

Table ^ lists the derived effective parameters. Only 
galaxies for which we also ha ve s pectroscopy are presented 
in this table. Tables Al 



and A 2 



contain the derived effec- 
tive parameters for all galaxies brighter than 19™ and within 
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Table 4. Photometric parameters 

NOT HST 

2 2 

Galaxy Field logrc <M>o ■n~~ <M>o <M>c {V - Ic) logrc </i>c i^— <A»>c </^>c {V - Ic) 

[arcsec] (r,) (V,) [arcsec] (/e) (r.) (Vi) 



A665-1150" 


1 


1 


012 


21 


82 


2 


9 


22 


67 


22 


96 


1 


54 


1 


056 


21 


97 


5 


2 


22 


82 


23 


11 


1.53 


A665-1168 


1 





333 


19 


78 


2 


1 


20 


63 


20 


92 


1 


54 





376 


19 


93 


1 


3 


20 


78 


21 


06 


1.52 


A665-1196 


1 


-0 


352 


17 


73 


9 





18 


58 


18 


86 


1 


51 


-0 


149 


18 


58 


7 


9 


19 


43 


19 


68 


1.47 


A665-2105'' 


2 





876 


22 


09 


29 


7 


22 


92 


23 


07 


1 


20 


1 


003 


22 


42 


18 


5 


23 


23 


23 


41 


1.24 


A665-2108 


2 


-0 


205 


19 


48 


1 


5 


20 


33 


20 


61 


1 


53 
























A665-2111 


2 


-0 


170 


18 


56 


6 


7 


19 


41 


19 


74 


1 


62 


-0 


Oil 


19 


17 


4 


1 


20 


02 


20 


29 


1.51 


A665-2144 


2 


-0 


189 


18 


83 


1 


6 


19 


68 


19 


96 


1 


52 


-0 


096 


19 


22 


1 


6 


20 


06 


20 


29 


1.41 


A2218-L118 


1 





051 


19 


20 


5 


5 


20 


06 


20 


39 


1 


62 





097 


19 


30 


3 


1 


20 


16 


20 


48 




A2218-L148 


1 


-0 


075 


18 


53 


3 


1 


19 


39 


19 


70 


1 


59 





067 


19 


03 


6 


4 


19 


89 


20 


20 




A2218-L198 


1 


-0 


061 


19 


14 


2 





19 


99 


20 


29 


1 


55 


-0 


143 


18 


81 


1 


9 


19 


67 


19 


96 




A2218-L244 


1 





934 


22 


01 


2 





22 


87 


23 


15 


1 


53 





930 


22 


02 


1 


4 


22 


88 


23 


16 




A2218-L341 


1 


-0 


124 


18 


27 


2 


6 


19 


13 


19 


42 


1 


56 


-0 


116 


18 


26 


4 


1 


19 


11 


19 


41 




A2218-L391" 


1 


1 


403 


23 


07 


3 


6 


23 


92 


24 


18 


1 


48 


1 


412 


23 


25 


21 





24 


10 


24 


37 




A2218-L430 


2 


-0 


033 


19 


06 


1 


4 


19 


92 


20 


24 


1 


62 


-0 


061 


18 


95 


1 





19 


81 


20 


13 




A2218-L482 


2 





036 


18 


77 


2 


2 


19 


63 


19 


95 


1 


62 
























A2218-L535 


2 





238 


19 


63 


1 


9 


20 


49 


20 


82 


1 


63 

























Notes. " Brightest cluster galaxy. ^ E+A galaxy. Identifications of galaxies in A2218 are from Le Borgne et al. (1992). Field - the NOT field 
that contains the galaxye, see Table Q for seeing values. The observed magnitudes labeled Ic and V and the colors {V — Ic) are standard 
calibrated to the Johnson-Kron-Cousins photometric system, and corrected for galactic extinction. The formal fitting uncertainties are 
typically as follows, logrc: ±0.006 for NOT data, ±0.002 for HST data. <At>e: ±0.02 for NOT data, ±0.007 for HST data. The typical 
(random) uncertainty on the {V — Ic) is 0.01. The goodness-of-fit is given as X'^/^pixi where Npi^ is the number of pixels within the 
fitting radius, = X/^*^' ~ '^i,model)^/'''|i with the sum over all the pixels within the fitting radius, rii is the signal in pixel i, n.i_niodol 
is the model value in the pixel, and tr? the noise (photon- and read-out-noise) in the pixel. 

Magnitudes are fully corrected Gunn r magnitudes in the rest frames of the clusters, see equations ^ and (^). Magnitudes Vz are 
fully corrected Johnson V magnitudes in the rest frames of the clusters, see equations M) and (ti). 



the fields covered by the NOT images. The magnitudes were 
standard cahbrated as described in Appendix A. The Ic 
magnitudes from the HST data have been offset to consis- 
tency with the ground-based magnitudes, see Appendix A. 
Table ^ gives typical values of the formal fitting uncertain- 
ties. However, a comparison of the NOT and the HST data 
results in more realistic estimates of the uncertainties, see 
Appendix A. We find the uncertainties of log re and </i>(, 
to be ±0.078 and ±0.29, respectively. Because of the cor- 
relation of the errors in logrc and <jj,>^, the combination 
logrc — 0.328</i>c that enters the FP has an uncertainty 
of only ±0.018. In the analysis of the FP for the clusters 
(Section 5), we use the average of the effective parameters 
determined from the HST observations and the NOT obser- 
vations, whenever results from both telescopes are available. 



2.2 Colors 

The colors of the galaxies were determined from the ground- 
based observations. First, the I-band images were convolved 
to the seeing of the V-band images. This minimizes the er- 
rors in the derived colors due to mismatched seeing in the 
two passbands. The colors were derived within apertures 
with radii l'.'64 (10 pixels) and 2'.'94 (18 pixels) for A665 
and A2218, respectively. The larger aperture for A2218 was 
necessary due to the very poor seeing of the V-band images. 
The colors within these apertures are in the following used 
as global colors of the galaxies, see Tab le H for the sample 
that also has spectroscopy and Tables Al and A2 for all 
galaxies brighter than 19™. For nearby E and SO galaxies. 



ally of the order 0^05 per dex in the radius, for colors based 
on Johnson B and Gunn r (e.g., j0rgensen et al. 1995a). 
The color gradients in the aperture magnitudes are smaller, 
typically less than 0^04 per dex in the radius. Assuming 
that the color gradients of the galaxies at redshift 0.18 are 
similar to those of nearby galaxies, the derived colors are ex- 
pected to deviate less than ±0.02 from global colors within 
the effective radii. This results in a negligible effect of less 
than 0^002 on the calibrated mag nitu de in Gunn r in the 



2.3 



rest frame of the clusters, cf. Sect. 

For the galaxies in A665 colors were also derived from 
the HST images. We used apertures with radii of l'.'5. The 
aperture size is smaller than those used for the NOT images 
in order to take advantage of the better spatial resolution 
of the HST images and limit the amount of contamination 
from neighboring objects. No correction for zero point off- 
sets between the ground-based photometry and the HST 
photometry were applied to the colors, since we have no 
reliable way of checking the zero point for the HST V mag- 
nitudes. Thus, we implicitly assume that the zero point dif- 
ference between the ground-based photometry and the HST 
photometry is the same in the I-band and the V-band. The 
colors derived from the HST observations are on average 
0.04 mag bluer than the ground-based colors, cf. Table m 
This may be caused by errors in the HST colors due to the 
uncertainties in the magnitude zero points. Alternatively, 
uncertainties may be introduced by the poor seeing of the 
ground-based V images. The effect on the magnitudes cali- 
brated to Gunn r in the rest frame of the cluster will be less 
than 0™003, cf. Sect. El . 



the color gradients in the local surface brightnesses are usu- 
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Table 5. Uncertainties affecting the photometric calibration 



Cluster 


Tel. 


Calibration/Parameter 


Random unc. Systematic unc. 




NOT 


Instrumental Jc 




0.01 


A 991 S 


MOT 


Instrumental Ic 




0.015 


A fif,'^ A 991 S 


NOT 


Standard Ic 


0.021 




A fkfi^ A 991 S 


NOT 


Standard V 


0.014 




A fkfi^ A 991 S 


NOT 


Standard {V - Ic) 


0.025 




A f\fi^ 

I\DOO 


HST 


Holtzman et al. (1995) Ic 


U.UIO 




A665 


HST 


Holtzman et al. (1995) {V - Ic) 


0.028 




A665 


HST 


Offset to ground-based Ic 




0.009 


A2218 


HST 


Holtzman et al. (1995) Ic 


0.017" 




A2218 


HST 


Offset to ground-based Ic 




0.013 


A665,A2218 


NOT,HST {V - /c)meas VS. {V - /c)r„ 


0.02 




A665,A2218 


NOT 


Total (V - Ic) 


0.032 




A665 


HST 


Total {V - Ic) 


0.034 




A665,A2218 


NOT, HST Gunn r in rest frame, transformation 


0.02'' 




A665,A2218 


NOT, HST Johnson V in rest frame, transformation 


0.02'' 




A665 


NOT 


Total 


0.028 


0.01 


A665 


HST 


Total 


0.024 


0.02 


A665 


NOT 


Total 


0.032 


0.01 


A665 


HST 


Total 


0.029 


0.02 


A2218 


NOT 


Total 


0.028 


0.015 


A2218 


HST 


Total 


0.024 


0.028 


A2218 


NOT 


Total 


0.032 


0.015 


A2218 


HST 


Total 


0.029 


0.028 



Notes. " Contributions from the transformation (V — Rc) = 0.52(y — Ic) and the transformation from 
Holtzman et al. (1995). '' Contribution from transformation, only. 



2.3 Calibration to the rest frame of the clusters 

In order to compare the FP for A665 and A2218 with results 
for nearby clusters we calibrate the photometry to Gunn r 
in the rest frame of the clusters. Following Frei & Gunn 
(1994), the magnitude in a broad passband can be written as 
rrii — Ci — 2.5 log F(i/i), with Ci = Abi —48.60. Abi axe given 
by Frei & Gunn. F{ui) is the flux density at the effective 
wavelength of the passband. We assume the flux density in 
the redshifted Gunn r band is related to the observed V- 
band and 7c-band as F{i'r{z)) = F{vv)°' F{vi^Y~°' , cf. van 
Dokkum & Franx (1996). The magnitudes in Gunn r in the 
rest frame of the clusters can then be written 

= Jc + a(V-/c)-ci^ -Q(cv-ciJ + Cr + 2.51og(l + z)(l) 

where a is derived from the spectral energy distribution. 
The Ic magnitudes and the colors {V — Ic) are the observed 
standard calibrated values, corrected for galactic extinction. 

A665 has z=0.181 (Oegerle et al. 1991). From the spec- 
tral energy distributions for E and Sa galaxies given by Cole- 
man, Wu & Weedman (1982) we find qa665 = 0.053. The 
resulting calibration to reach Gunn r in the rest frame is 

= ic + 0.053(V - Ic) + 0.77 (2) 

A2218 has z=0.176 (Le Borgne, Pello & Sanahuja 1992). We 
find aA22i8 = 0.062, and the resulting calibration is 

= Jc + 0.062(1/ - Ic) + 0.76 (3) 

The mean surface brightnesses, </i>g, calibrated to Gunn r 
in the rest frame of the clusters are given in Table ^. 

Table ^ also lists the mean surface brightnesses cali- 
brated to Johnson V in the rest frame of the clusters. The 
calibrations to reach Jonhson V are 



V,=Ic + 0.48(V - Ic) + 0.40 (4) 
for A665 and 

= Je + 0.49(1/ - Ic) + 0.39 (5) 
for A2218. 

Because the photometry for the best comparison sample 
(the Coma cluster) is available in Gunn r, and also because 
of the uncertainty in the derived colors, cf. Section 2.2, we 
will primarily use the mean surface brightnesses calibrated 
to Gunn r throughout the rest of this paper. 

2.4 Summary of uncertainties 

Table |^ summarizes the different sources on uncertainty af- 
fecting the photometric calibration. The entries noted as 
"Total" uncertainty give the total uncertainty on the listed 
parameter, while all other entries give the contribution from 
the individual calibration step. The uncertainties on the in- 
strumental Ic magnitudes for the NOT data refer to the 
consistency of the magnitude zero points of the co-added 
images, cf. Sect. Al. The uncertainties on the standard cal- 
ibrated magnitudes for the NOT data are the rms scatter 
of the standard transformations, cf. Sect. A3. For the HST 
data the uncertainties consist of random uncertainties due to 
the uncertainties in the transformations adopted from Holtz- 
man et al. (1995), cf. Sect. A3, and systematic uncertainties 
in matching the HST photometry to the ground-based sys- 
tem. The colors {V — Ic) are subject to an additional uncer- 
tainty because we measure the colors within fixed aperture 
sizes rather than within the effective radii, cf. Sect. 2.2. Fi- 
nally, the transformations used to calibrate the magnitudes 
to Gunn r and Johnson V in the rest frames of the clus- 
ters contribute to the uncertainties. The last part of Table 
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^summarizes the total uncertainties, random or systematic, 
on the cahbrated magnitudes and Vz- The magnitudes 
have typical random uncertainties of 0™03, and typical sys- 
tematic uncertainties of CP02. 



3 THE SPECTROSCOPY 

A665 galaxies were selected for the spectroscopic observa- 
tions from a ground-based R-band image kindly provided 
by M. Birkinshaw. The flux of the galaxies within a I'.'S x 8" 
aperture was used as selection criterion. Spectra of 7 galax- 
ies in A665 were obtained with the MMT in January 1991, 
cf. Table |^. The total integration time was 9 hours. For one 
of the galaxies the signal-to-noise was too low to derive a 
reliable value of the velocity dispersion. 

For A2218 we used the catalog by Le Borgne, Pello & 
Sanahuja (1992) as the main source for the selection of galax- 
ies, supplemented with the catalog by Butcher, Oemler & 
Wells (1983). Galaxies were selected on magnitude and color. 
The (B — R) color was required to be within 0.15 of the color- 
magnitude relation for the red cluster galaxies. One redder 
galaxy, BOW97, was included. This galaxy turned out to be 
a background galaxy. Spectra of 10 galaxies in A2218 were 
obtained with the KPNO 4m in June 1994, cf. Table |l|. The 
total integration time was 11.8 hours. 

The basic reductions of the spectroscopy are described 
in Appendix B, which also shows the final spectra of the 
galaxies. Here we describe the determination of the velocity 
dispersions (see also Franx 1993ab for the galaxies in A665) . 

3.1 Instrumental resolution and template stars 

Determination of velocity dispersions of galaxies is usually 
done with spectra of template stars obtained with the same 
instrumentation as the galaxy spectra. Due to the large red- 
shifts of A665 and A2218 this will not work for these clus- 
ters. A further complication is that the spectra of A665 and 
A2218 were obtained through multi-slit masks and the in- 
strumental resolution varies with wavelength and from slit- 
let to slit-let. 

Following van Dokkum & Franx (1996), we mapped the 
instrumental resolution from He-Ne- Ar lamp exposures. Fur- 
ther, the instrumental resolution was determined from the 
sky lines in the final average spectra of the galaxies. In all 
cases Gaussian profiles were fitted to the emission lines. Fig. 
|l| shows the instrumental resolution versus wavelength for 
representative slit-lets for the A2218 spectra. The determi- 
nations from lamp exposures and from sky lines agree. Low 
order polynomials were fitted to the instrumental resolution 
versus wavelength for each slit-let. 

We used a spectrum of HD192126 as the template star 
for the A665 spectra. The spectrum of the star was ob- 
tained with the Canada-France-Hawaii Telescope by K. Kui- 
jken. The wavelength range for the spectrum is 4650-5650A, 
with a resolution of crinst = 0.56A. As template stars for 
the A2218 spectra we used spectra of HD56224, HD72324 
and HD 102494 obtained with the Keck Telescope in May 
1996 by D. Kelson. The wavelength range of these spectra 
is 4002-5314A, with a resolution of ainst = 0.9A. 

For each galaxy spectrum the spectra of the template 
stars were convolved to the instrumental resolution of the 




4500 5000 5500 6000 6500 
Wavelength [A] 



Figure 1. The instrumental resolution for the A2218 spectra, 
as derived from He-Ne- Ar calibration spectra. Triangles - central 
slit-let; boxes and circles - the two slit-lets closest to the edges of 
the CCD. The filled symbols are based on sky lines. The resolution 
derived from the sky lines agree with the resolution derived from 
the calibration spectra. The solid lines show the fitted functions. 
The dot-dashed lines show the approximate wavelength interval, 
which was used for determination of the velocity dispersions. 

galaxy spectrum. The variation of the resolution with wave- 
length was taken into account. 

3.2 Fourier Fitting 

The velocity dispersions, a, were determined by fitting the 
galaxy spectra with the template spectra using the Fourier 
Fitting Method (Franx, lUingworth & Heckman 1989b). The 
Fourier transforms of the spectra were filtered to avoid the 
low frequencies. We used a limit equivalent to ~(100A)~^ 
in the rest frame of the clusters. This is the same limit used 
by j0rgensen, Franx & Kjaergaard (1995b). In the case of 
A2218, the final velocity dispersion is the average of the de- 
terminations with the three different template stars. The for- 
mal fitting uncertainty on log a is typically 0.023. For A665 
the mean spectra over the central l'.'8 were fitted. Thus, the 
aperture size is l'.'8xl'.'8. The aperture size used for A2218 
was l'.'9x3'.'45. 

For the A2218 galaxies we experimented with fitting 
the spectra in a limited wavelength range around the G- 
band region (rest frame 4080-4700A) and around the Mgb 
region (rest frame 4800-5310A), as well as the full available 
wavelength range (rest frame 4080-5310A). Fig. | shows the 
comparison of the velocity dispersions derived from the G- 
band region and from the Mgb region. The median offset in 
log a is 0.00 ± 0.02. In the following the results from fitting 
the full wavelength range are used. 

The velocity dispersions were aperture corrected to a 
circular aperture with diameter 1.19h~^kpc, equivalent to 
3.4 arcsec at the distance of the Coma cluster. We adopt 
the technique for the aperture correction established by 
j0rgensen et al. (1995b). The applied aperture corrections 
are 0.0215 and 0.028 in log a for A665 and A2218, respec- 
tively. Table ^ gives the raw determinations of the velocity 
dispersions as well as the aperture corrected values. 

The velocity dispersion measurements for A665 and 
A2218 are based on two independent observing runs. Fur- 
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Figure 2. The difference between velocity dispersions derived 
from the Mgb region (rest frame 4800-5310A) and from the G- 
band regon (rest frame 4080-4700A) versus the velocity disper- 
sions derived from the full wavelength range. Only data for galax- 
ies in A2218 are shown. The median offset in log c is 0.00 ± 0.02. 



ther, in the analysis (Sect. 5) we will use the data together 
with measurements from other independent observing runs. 
We have no way of testing the consistency of the velocity dis- 
persion measurements from the different runs, since there 
are no galaxies in common between the runs. For nearby 
galaxies (redshifts less than 0.1) it is in general found that 
velocity dispersion measurements from different runs can de- 
viate with 0.01 to 0.025 in logo- (e.g., JFK1995b; Smith ct 
al. 1997; Wegner et al. 1998). We adopt the mean offset of 
0.023 in log CT found by Wegner et al. in their analysis of 33 
observing runs to be the typical systematic uncertainty for 
velocity dispersion measurements from independent observ- 
ing runs. The uncertainty on the applied aperture correction 
is expected to be significantly smaller than this systematic 
uncertainty. 



Table 6. Spectroscopic parameters 



Galaxy log fx log cr crioo ,T ^ S/N 



A665-1150 


2.447 


2.468 


0.019 


0.1831 


54.5 


A665-1168 


2.356 


2.378 


0.019 


0.1785 


54.5 


A665-1196 


2.334 


2.355 


0.028 


0.1691 


32.6 


A665-2105 


2.127 


2.149 


0.035 


0.1699 


40.7 


A665-2108 








0.1939 


16.5 


A665-2111 


2.369 


2.391 


0.025 


0.1860 


37.2 


A665-2144 


2.267 


2.288 


0.034 


0.1729 


29.4 


A2218-L118 


2.280 


2.308 


0.023 


0.1748 


32.2 


A2218-L148 


2.371 


2.399 


0.023 


0.1833 


27.6 


A2218-L198 


2.309 


2.337 


0.027 


0.1737 


22.4 


A2218-L244 


2.287 


2.315 


0.020 


0.1772 


32.2 


A2218-L341 


2.451 


2.479 


0.023 


0.1630 


25.2 


A2218-L391 


2.483 


2.511 


0.022 


0.1731 


38.1 


A2218-L430 


2.093 


2.121 


0.035 


0.1810 


21.7 


A2218-L482 


2.342 


2.370 


0.018 


0.1788 


35.2 


A2218-L535 


2.270 


2.298 


0.021 


0.1809 


31.1 


A2218-BOW33 


2.186 


2.214 


0.038 


0.1734 


18.2 


A2218-BOW97 








0.291 


12.7 



Notes, (obs) — raw measurements. (3'.'4) - aperture corrected to a 
circular aperture with diameter 1.19 h^^ kpc, equivalent to 3'.'4 at 
the distance of the Goma cluster, ciog „ - formal uncertainty from 
the Fourier Fitting Method. S/N - signal-to-noisc per Angstrom 
in the wavelength interval used for the Fourier fitting. Identifica- 
tions for A2218 are from Le Borgne et al. (1992) and Butcher et 
al. (1983). 

the coefficients of the FP we have also used data for galax- 
ies in the Hydral cluster. The main purpose of including 
the Hydral cluster is to add more high luminosity galax- 
ies to the nearby sample. We use the sample of E and SO 
galaxies by Milvang- Jensen & j0rgensen (1999). The sam- 
ple covers the central 68' x 83' of the cluster. A total of 45 
galaxies have photometry (in Gunn r) and spectroscopy. The 
sample is 80% complete to niT = 14™ 5 in Gunn r (equiva- 
lent to Mr^ = -20™0, for a distance modulus of 34™ 5). For 
the galaxies brighter than Mr^ = -20^75, the FP for the 
Hydral cluster is the same as for the Coma cluster. 



4 DATA FOR OTHER CLUSTERS 

4.1 The Coma cluster and the Hydral cluster 

The Coma cluster is in the following used as the low redshift 
reference for our study of changes in the FP as a function of 
the redshift. We use the sample of E and SO galaxies, which 
has photometry in Gunn r from j0rgensen et al. (1995a), see 
also j0rgenscn & Franx (1994). The sample covers the cen- 
tral 64' X 70' of the cluster. Velocity dispersions published by 
Davies et al. (1987), Dressier (1987), and Lucey et al. (1991) 
have been calibrated to a consistent system and aperture 
corrected to a circular aperture with diameter 1.19 h~^ kpc, 
equivalent to 3.4 arcsec at the distance of the Coma clus- 
ter, see j0rgensen ct al. (1995b). Further, we use new ve- 
locity dispersion measurements from j0rgensen (1999). A 
total of 116 galax;ies have available photometry and spec- 
troscopy. The sample is 93% complete to a magnitude limit 
of rriT = 15™05 in Gunn r, equivalent to M^^, = -20™75 
(for a distance modulus of 35™8, Ho = 50 kms"'^ Mpc~^ and 
qo = 0.5). 

For the purpose of testing for possible differences in 



4.2 Intermediate redshift clusters 

In the following analysis wo also use data for the three inter- 
mediate redshift clusters CL0024-I-16 (van Dokkum & Franx 
1996), CL1358-I-62 and MS2053-04 (Kelson et al. 1997). Van 
Dokkum & Franx and Kelson et al., who previously stud- 
ied these clusters, give effective radii, magnitudes Vz cal- 
ibrated to the V-band in the rest frame of the clusters, 
colors, and velocity dispersions. We establish the transfor- 
mations from the observed Ic magnitudes and (Re — Ic) to 
Gunn g in the rest frames of the clusters, using the tech- 
nique described in Section 2.3. We then use the relation 
{V -g) = -0.517(5 - r) - 0.027 (J0rgensen 1994), based on 
standard stars, to transform (14 — flz) to [g^ — r^), where 
Qz and Tz are Gunn g and Gunn r, respectively, in the rest 
frames of the clusters. This results in calibrations that allow 
us to transform from the published magnitudes Vz and the 
color (Rc — Ic) to Gunn r in the rest frames of the clusters. 
For CL0024-I-16 at a redshift of 0.39 we find 

- rz = 0.58(i?c - h) - 0.21 (6) 
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Figure 3. The FP edge-on for Coma, A2218, A665, CL1358-I-62, CL0024-I-16, and MS2053-04. The sources of the data are given on the 
panels. The skeletal symbols on panel (c) and (d) are the E+A galaxies. The photometry is calibrated to Gunn r in the rest frames of 
the clusters. The mean surface brightness log </>c = —OA{<fi>^ — 26.4) is in units of Lq/pc'^. The photometry is not corrected for the 
dimming due to the expansion of the Universe. The effective radii are in kpc {Hq = 50 kms~^ Mpc~^ and Qo = 0.5). The solid lines are 
the FPs with coefficients adopted from JFK96 and zero points derived from the data presented in the figure. Typical error bars are given 
on the panels; the thin and thick error bars show the random and systematic uncertainties, respectively. 



For CL1358+62 at z = 0.33 we find 

K -r^ = 0.50(i?c -/c) - 0.19 (7) 

The data for CL1358+62 include two E+A galaxies (Kelson 
et al.). For MS2053-04 at z = 0.58 we find 

K -r, =0.30(i?c -/c) "0.15 (8) 

The uncertainties on derived using these transforma- 
tions consist of the uncertainty due to the transformation 
between {V — g) and (g — r) (0.015) and the uncertainty due 
to the transformation to g^, (0^02), a total of 0'"025. We do 
not include additional uncertainty from the authors' orig- 
inal transformation to Vz since that uncertainty is closely 
correlated with the uncertainty on the transformation to g^. 
The uncertainty of 0™025 does not include the uncertainties 
in the standard calibrations to Ic and Rc established by van 
Dokkum & Franx and Kelson et al. 

The typical values of Vz — derived from equations (^ 
to ^ are 0.2. This is comparable to the average {V — r) color 
of nearby E and SO galaxies (e.g, j0rgensen et al. 1995a). 
Others found similar small color evolution for E galaxies in 
clusters with redshifts less than 0.6, e.g., Rakos & Schombert 



(1995) who used passbands shifted to match the redshift of 
the studied clusters. 



5 THE FUNDAMENTAL PLANE 

Fig. H shows the FP edge-on for A665 and A2218 as logcr 
versus log To + 0.82 log <7>e. The mean surface brightness 
log<7>c = -OA{<ii>^ - 26.4) is in units of Lo/pc^. The 
figure also shows the FP for the Coma cluster and for the 
three intermediate redshift clusters CL1358-I-62 [z = 0.33), 
CL0024-hl6 {z = 0.39) and MS2053-04 {z = 0.58). The 
E-l-A galaxies in A665 and CL1358-I-62 are excluded from 
the following analysis, but are shown on the figures. 

We adopt the coefficients for the FP from JFK96, 
(q,/3) = (1.24 ± 0.07,-0.82 ± 0.02). The coeflicients were de- 
rived for a sample of 226 galaxies in 10 nearby clusters. 
Photometry in Gunn r was used. The relations shown on 
Fig. ^ are FPs with these coefficients. 

The error bars on Fig. ^ showing the systematic errors 
include the systematic uncertainties in the photometry (cf. 
Table ^ and a 2% systematic uncertainty in re</>e° *^ due 
to any deviations from r^^* profiles and possible differences 
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Figure 4. The FP for Coma, Hydral and the five intermediate redshift clusters. Small triangles — galaxies in Hydral; other symbols as 
on Fig. ^ The photometry is calibrated to Gunn r in the rest frames of the clusters and corrected for the dimming due to the expansion 
of the Universe. The effective radii are in kpc. Hg = 50 kms"'^ Mpc"'^ and go = 0.5 were used, (a) The FP face-on. The dashed line 
marks the luminosity limit for the completeness of the Coma cluster sample, Mrrp = — 20^"75. The dotted line is the so-called exclusion 
zone for nearby galaxies, y —0.54a; -|- 4.2, first noted by Bender et al. (1992). (b) and (c) The FP edge-on. The solid line is the FP for 
the Coma cluster with coefficients adopted from JFK96. The small error bars on panels (a) and (b) refer to Coma and Hydral, while 
the larger error bars refer to the intermediate redshift clusters. The error bars on panel (c) refer to Coma and Hydral. The error bars 
for the intermediate redshift clusters cannot be summarized in a simple fashion, and we refer to Fig. M for their sizes. 



originating from differences in the fitting methods for the 
nearby and intermediate redshift galaxies (cf. Kelson et al. 
1997). The systematic uncertainties on log cr is shown as 
0.023 for the intermediate redshift clusters and zero for the 
Coma cluster, since the relevant uncertainty is the possi- 
ble inconsistency between the velocity dispersions for Coma 
galaxies and the intermediate redshift galaxies. 

Fig. ^ shows the FP face-on and in two edge-on views. 
The Hydral sample is included on this figure. The mean 
surface brightnesses have been corrected for the dimming 
due to the expansion of the Universe. The effective radii are 
in kpc (Ho = 50 kms~^ Mpc~^ and go ~ 0.5). 

The scatter around the FP is similar for A665 and 
A2218 (0.060 and 0.112 in log re) and for the Coma clus- 
ter (0.092 in logTe). Using the photometry calibrated to 
Gunn r, we find the scatter for CL1358-H62, CL0024-f 16 
and MS2053-04 to be 0.066, 0.064 and 0.076, respectively 
(see also van Dokkum & Franx 1996; Kelson et al. 1997). 
This is not significantly different from the scatter for Coma, 
A665 and A2218. The scatter for the Coma cluster reported 
here is somewhat larger than previous estimates based on 



smaller samples (e.g., j0rgensen et al. 1993, 1996; Lucey et 
al. 1991). 

5.1 The coefficients of the FP 

The sample of intermediate redshift galaxies is heavily bi- 
ased towards intrinsically bright galaxies. The faintest galax- 
ies observed in the intermediate reshift clusters typically 
have Mi-rp = — 21™65. The selection effect is clearly visible 
in the face-on view of the FP, Fig. ^a. Further, the edge-on 
view on Fig. ^ shows that the samples are limited in both 
the velocity dispersions and in logro -I- 0.82 log <I>c. The 
selection effects make it difficult to test for possible differ- 
ences between the coefficients of the FP for Coma and the 
FP for the intermediate redshift clusters. 

Fig. ^ shows the FP as the relation between the M/L 
ratio and the mass. The masses were derived as M = 
Sr^a^/G (cf.. Bender, Burstein & Faber 1992). This fig- 
ure and Fig. Wc indicate that the coefficient a for the log a 
term in the FP is different for the intermediate redshift clus- 
ters and for the Coma cluster. This was also noted by van 
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Dokkum & Franx (1996) in their analysis of the CL0024+16 
observations. 

We tried to fit the FP to the intermediate redshift clus- 
ters individually, in order to test if the coefficients of the FP 
for the intermediate redshift clusters are different from those 
for nearby clusters. A total of 35 galaxies were included in 
the fits. The E+A galaxies in A665 and CL1358-I-62 were 
omitted, since we aim at fitting the FP for the E and SO 
galaxies, only. The fits were determined by minimization 
of the sum of the absolute residuals perpendicular to the 
relation. The uncertainties were derived using a bootstrap 
method (cf., JFK96). All the fits gave /3 values in the inter- 
val —0.75 to —0.81. The a values are in the interval 0.54- 
1.74 and, due to the small sample sizes, the uncertainties 
are large. We therefore repeated the fits with f3 = —0.79 
for all the clusters. This gives a in the interval 0.45-1.42, 
with uncertainties of 0.14-0.70. No significant change of a 
with redshift was detected. It must await larger samples of 
galaxies to map the possible change of the coefficients as a 
function of the redshift. 

Next, we fit the FP to the five clusters as parallel 
planes, under the assumption that the FPs for these clusters 
have the same coefficients. For the photometry calibrated to 
Gunn r, this gives (a,/3) = (0.89±0. 17,-0. 79±0.04). In John- 
son V we find (a,/3)=(0.92 ± 0.13,-0.79 ± 0.04). A fit to the 
Coma cluster galaxies alone gives (ci,/3) = (1.27±0.08,— 0.81± 
0.03), for photometry in Gunn r. The difference between a 
for the Coma cluster sample and a for the intermediate red- 
shift sample is 0.38. Of the 1000 bootstrap fits derived for 
the FP for the intermediate redshift sample, 3.8 percent have 
a values that deviate from a — 0.89 with more than 0.38. 
Hence, the difference in a is significant at the 96 percent 
level. 

In order to explore the effect of the different selection 
criteria we repeated the fit to the Coma cluster galaxies first 
enforcing a magnitude limit of M,.rj, = —21™ 65, and next 
enforcing limits of log ct > 2.12 and log ro -I- 0.82 log <7>o > 
2.80. Applying the magnitude limit, does not give a result 
significantly different from the fit for the full sample, we find 
(a,/3)=(1.30 ± 0.17,-0.82 ± 0.06). The limits in logo" and 
logrc + 0.82 log </>c resuh in (q,/3)=(1.09 ± 0.27,-0.85 ± 
0.05). Though this is formally not different from the fit to the 
full sample, the smaller value of the a coefficient indicates 
that part of the difference in a between the Coma cluster 
sample and the sample in the intermediate redshift clusters 
may be due to selection effects. 

To further test if the difference in a may be due to selec- 
tion effects, we have performed two simulations that are both 
based on fitting the FP to random sub-samples of galaxies 
selected from the Coma cluster sample plus the Hydral clus- 
ter sample. The sub-samples each contain 35 galaxies in or- 
der to match the size of the intermediate redshift sample. In 
the first simulation, the sub-samples were selected such that 
their luminosity distributions match the luminosity distri- 
bution for intermediate redshift cluster sample. In the sec- 
ond simulation, the sub-samples were selected to match the 
distribution of logr-o + 0.82 log </>o. The mean luminosity 
evolution of the galaxies in the intermediate redshift clusters 
(see Section ^.2| ) was taken into account in both simulations. 
However, the results do not depend on whether this correc- 
tion is applied. 

We have included the Hydral cluster sample in these 
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Figure 5. (a) The M/L ratio as function of the mass. The dashed 
line marks the luminosity limit for the completeness of the Coma 
cluster sample, M^^ = -20™75. (b) The M/L ratio as function 
of the velocity dispersion. Solid line — relation fitted to the Coma 
cluster; dotted-dashed line - relation fitted to the intermediate 
redshift clusters. Symbols as in Fig. |^. The small and large error 
bars refer to the random uncertainties aff'ecting the Coma cluster 
and the intermediate redshift clusters, respectively. 



simulations because the Coma cluster sample has only six 
galaxies brighter than Mr^p = —23™. Our sample of galax- 
ies in the intermediate redshift clusters include 13 galaxies 
brighter than Mr.p = —23™. Thus, using only the Coma clus- 
ter sample for the test could potentially bias the result, since 
the bright galaxies would appear several times in each sub- 
sample. The combined Coma and Hydral sample contains 
10 galaxies brighter than A-I,-^ = —23™. While this may not 
be enough to eliminate all bias due to some of the bright 
galaxies being included more than once in each sub-sample, 
we expect the bias to be fairly small. 

Figs, ^a and c show the distribution of a values for 1000 
random sub-samples for each of the two simulations. The 
values for the full Coma cluster sample and for the inter- 
mediate redshift sample are overplotted. Eight (0.8 percent) 
of the sub-samples selected to match the luminosity distri- 
bution result in values of a within the one sigma error bar 
of the coefficient found for the intermediate redshift sample. 
For the sub-samples selected to match the distribution in 
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Figure 6. (a) Histogram of tlie values of tlie coefficient a of the FP for 1000 random sub-samples of galaxies selected from the Coma 
cluster and Hydral cluster sample. The sub-samples are selected so their luminosity distributions match the luminosity distribution for 
the intermediate redshift sample. The values found for the intermediate redshift sample and for the full Coma cluster sample are also 
shown, (b) Histogram of the difference between the a values derived by bootstrapping the full Coma sample and the a values derived from 
the random sub-samples selected to match the luminosity distribution. The absolute difference |Aa| is larger than 0.38 for 1.8 percent 
of the cases, (c) Histogram of the values of the coefficient a of the FP for 1000 random sub-samples of galaxies selected from the Coma 
cluster and Hydral cluster sample. The sub-samples are selected so their distribution in log re -|- 0.82 log </>e match the distribution for 
the intermediate redshift sample. The values found for the intermediate redshift sample and for the full Coma cluster sample are also 
shown, (b) Histogram of the difference between the a values derived by bootstrapping the full Coma sample and the a values derived 
from the random sub-samples selected to match the distribution in logro + 0.82 log </>c. The absolute difference |Ao| is larger than 
0.38 for 3.8 percent of the cases. 



log Te -I- 0.82 log </>c, 20.3 percent give a values within the 
one sigma error bar. 

The distributions of the /3 values for the sub-samples en- 
compass for both simulations the value for the Coma cluster 
sample as well as the value for the intermediate redshift sam- 
ple, which in turn are not significantly different from each 
other. 

In order to construct two-sided tests, we compare the a 
values derived from the two simulations of the sub-samples 
to those derived by bootstrapping the full Coma sample. 
The comparison with the first simulation gives a measure 
of how often we can expect two samples with luminosity 
distributions similar to the intermediate redshift sample and 
the Coma cluster sample, respectively, and drawn from the 
same underlying distribution, to give a difference in a larger 
than the measured difference of 0.38. The comparison with 
the second simulation gives the same measure for samples 
selected on the distribution in logrc 4- 0.82 log </>o. Fig. 



and d show the distributions of the a difference for the 
two comparisons. For the sub-samples selected to match the 
luminosity distribution, 18 of the 1000 cases, or 1.8 percent, 
have an absolute difference |Aq| larger than 0.38. For the 
sub-samples selected to match the distribution in logro + 
0.82 log </>o, 3.8 percent of the cases have |Aci| > 0.38. 

Based on these simulations, we conclude that the differ- 
ence in the a coefficient of the FP is not likely to be due to 
the difference in the luminosity distributions of the two sam- 
ples. It is possible that part of the difference in the a coeffi- 
cient is due to the distribution of log re+0.82 log <7>e being 
different for the two samples. The different appearances of 
the distributions shown on Figs. |^ and c add to the evidence 
that the selection effect present in log r-o -f 0.82 log </>c may 
affect the derived value of a. 

The coefficients we find for the FP for the intermediate 
redshift clusters imply that 

tr/T 0.27 0.87 »,f0.44 -0.17 /n^ 

M/L (X ro a oc M (9) 



|3 
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Table 7. Uncertainties on the FP zero point 



■ Clusters ■ 



Source 




Coma 


A665 


A2218 


random / systematic 


Photometric zero point cr{<fj,>^) 




0.015 


0.02 


0.02 


systematic 


Velocity dispersion 




0.000 


0.023 


0.023 


systematic 


Coma peculiar velocity, 300km/s, 


o-(logro) 


0.018 






systematic 


Total systematic uncertainty, (t(1o 




0.023 


0.035 


0.035 


systematic 


Uncertainty on the FP zero point' 


(T(logre) 


0.009 


0.027 


0.037 


random 



Note. The uncertainty on the FP zero point is crgt/Af^/^ where trgt is the rms scatter of the FP and 
A'^ is the number of observed galaxies. 



This should be compared with M/L oc 
^0.24^^-0.02 j-Q^jj^ nearby clusters (JFK96). The diflter- 
ence may indicate that the low mass galaxies show a stronger 
luminosity evolution than the more massive galaxies. The 
driving factor may however be the velocity dispersion rather 
than the mass of the galaxies. A direct fit to the M/L ratio 
as a function of the velocity dispersion gave for the nearby 
clusters M/L oc ct" **^ (JFK96). Fig. |b shows the M/L ra- 
tio versus the velocity dispersion for the Coma cluster and 
the intermediate redshift clusters. The different slopes for 
the two samples is clearly seen. For the Coma cluster alone 
we find M/L oc q-0-(>6±o.i3 ^ while the sample of galaxies in 
the intermediate redshift clusters gives M/L oc cri-49±0-29 
The relations for the Coma cluster and for the intermediate 
redshift clusters are overplotted on Fig. |^. 

If real, the difference in the slopes of the FP (and the 
slopes of the M/L-ct relation) may be due to bursts of star 
formation in the low mass galaxies in the intermediate red- 
shift clusters leading to a lower M/L ratio for these galaxies. 
Alternatively, the samples in the intermediate redshift clus- 
ters may contain some early-type spirals. 



5.2 The evolution of the M/L ratio 

The zero point of the FP depends on the cosmological eff'ects 
(surface brightness dimming, the value of go, and the value 
of the cosmological constant), and the evolution of the galax- 
ies. We assume a cosmological constant of zero. We adopt 
the coefficients for the FP for nearby clusters (JFK96) for 
all the clusters. Because of the possible difference in the co- 
efficients of the FP (cf. Section 5.1) the following discussion 
of the evolution of the M/L ratios should be understood as 
the average evolution of the available samples of galaxies. 
Specifically, the results refer to galaxies with masses above 
6.5 • 1O^°M0 {Ho = 50 kms"^ Mpc"^ and go = 0.5). 

Table M summarizes the sources of uncertainty afi'ecting 
the zero points of the FP for Coma, A665 and A2218. The 
total systematic uncertainty on the FP zero point is the sum 
of the systematic uncertainties and is given as a{\ogrc). The 
random uncertainties introduced by the random photomet- 
ric uncertainties and fitting uncertainties are included in the 
uncertainty on the FP zero point since this is based on the 
measured rms scatter of the FP. 

On the following figures the error bars reflect the uncer- 
tainties listed in Table |^. The uncertainties for the other in- 
termediate redshift clusters were adopted from van Dokkum 
& Franx (1996) and Kelson et al. (1997), and where neces- 
sary the total uncertainties were derived in the same way as 
for A665 and A2218. 



V 




0.05 0.1 0.15 0.2 0.25 
log (1 + z) 



Figure 7. The offsets in the surface brightnesses derived from 
the FP versus the redshift of the clusters, qo = 0.5 was used. 
The small and large error bars show the random and system- 
atic uncertainties, respectively^These were derived as (T(</i>g) = 
(7(logrc) • 2.5/0.82. See Table | for values of o-(logrc). The solid 
line is the expected relation for an expanding Universe and no 
evolution of the galaxies. The data are in agreement with an ex- 
panding Universe and a slow luminosity evolution. A smaller qo 
results in smaller values of A<fi>, and therefore a larger lumi- 
nosity evolution. 



Fig. ^ shows the offsets in the surface brightnesses, 
A<^>, based on the zero points of the FPs. We used 
go = 0.5. No correction for the expansion of the Universe 
has been applied and the solid line on the figure shows the 
expected relation for an expanding Universe. A smaller go 
results in smaller values of A</i>. The data are consistent 
with an expanding Universe and a slow luminosity evolution 
of the galaxies. In the following we correct the data for the 
expansion of the Universe, and then use the FP to study the 
mean evolution of the galaxies in the clusters. 

The FP zero point difi^erences between the intermediate 
redshift clusters and the Coma cluster can be interpreted as 
differences in the M/L ratios. With the FP coefficients from 
JFK96 the M/L ratios can be expressed as 

log M /L = 0.49 log a + 0.22 log - 7/O.82 + est 

0.02 log ro 



0.49 log CT + 0.22 log ro 
0.24(2 log a -Flog To) - 
-7/0.82 + est 



(10) 



(cf., JFK96), where 7 is the FP zero point. Under the as- 
sumption that no merging or accretion takes place, the mass 
and therefore (2 log a + log r^) is constant. We can then de- 
rive the differences in the M/L ratios from the differences in 
the FP zero points, ignoring the small term 0.02 log Tc. Fig. 
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Figure 8. The evolution of the M/L ratio as function of the 
redshift of the clusters. The photometry has been calibrated to 
Gunn r in the rest frame of the clusters. The small and large error 
bars show the random and systematic uncertainties, respectively^ 
These were derived as a(\og M / L) = o-(log ro)/0.82. See Table j?] 
for values of (T(logre). (a) go = 0.5; (b) qo = 0.15. The top axis 
on each panel shows the time of emission in units of the present 
age of the Universe. Solid lines - least squares fits to the data. 
Dashed lines - passive evolution with Zform = °Oi 5 S'ld 1 (top to 
bottom). 



|8| shows the differences in the M/L ratios as a function of 
the redshifts. 

The galaxies in A665 and A2218 have on average M/L 
ratios (11 ± 7)% smaller than the M/L ratios of the galaxies 
in the Coma cluster. Together with the previous data from 
van Dokkum & Franx (1996) and Kelson et al. (1997), this 
shows the gradual and slow evolution of the M/L ratio of 
the galaxies as a function of the redshift. A least squares fit 
to the all data shown on Fig. ^ gives A log MjL,. = (— 0.26± 
0.06) Az for (?o = 0.5 and AlogAf/L^ = (-0.34 ± 0.06)A2 
for f/o = 0.15. For the photometry calibrated to Johnson V 
we find A log M/L^/ = (-0.29±0.08)Az and AlogM/L^ = 
(-0.37 ± 0.08)A2 for go = 0.5 and 0.15, respectively. Our 
results in Johnson V are in agreement with the results by 
van Dokkum & Franx and Kelson et al. 

The M/L ratio of a stellar population can be modeled 

as 



MjL oc 



(11) 



where t is the age of the stellar population (Tinsley 1980). 
The coefficient k depends on the initial mass function (IMF) 
of the stars and on the passband. Based on the single burst 
stellar population models from Vazdekis et al. (1996), we 
find for photometry in Gunn r that ~ 0.9 — 0.2a; where x 



is the slope of the IMF. A Salpeter (1955) IMF has x = 1.35 
(equivalent to K,r ~ 0.6). For photometry in Johnson V, 
the expected evolution is nearly the same. The models from 
Vazdekis et al. give nv ~ 0.95 — 0.2a;. 

Equation ( |ll[ ) can be used to relate the change in the 
M/L ratio with redshift to the formation redshift, Zform, and 
the value of go- The dashed lines on Fig. |^ show the ex- 
pected relations for k = 0.6 and Zform = oo,5 and 1. Models 
with shallower IMFs make the slope of the expected relation 
steeper, while a model with steeper IMFs make the slope of 
the relation shallower. 

For go = 0.5, the models with a Salpeter IMF and 
Zform > 5 are consistent with the data within 2g. A for- 
mation redshift of Zform = 2.4 represents a 3(t deviation. 
Models with Zform smaller than 2.4 deviate more than 3(t 
from the fit to the data. If we assume go = 0.15 more mod- 
els are consistent with the data. A model with a Salpeter 
IMF then implies Zform ~ 2.5. Models with Zform in the in- 
terval 1.7-6.5 are consistent with the data within 2a. Models 
with Zform < 1.4 deviate more than 3(t from the fit to the 
data. 

Because of the larger number of clusters analyzed here 
we are able to put stronger formal limits on the formation 
redshift, zform, than done in previous studies (van Dokkum 
& Franx 1996; Kelson et al. 1997). However, the derived 
model constraints should only be taken as rough guidelines. 
The correct interpretation of the data is most likely rather 
more complicated than indicated here. The correct value of 
K is not known. Further, the evolution of E and SO galax- 
ies cannot be viewed as a single burst event. The presence 
of younger stellar populations in the galaxies would imply 
stricter lower limits on Zform- The most fundamental assump- 
tion for the interpretation of the data is that the observed 
E and SO galaxies in the intermediate redshift clusters in 
fact evolve into galaxies similar to the present day E and SO 
galaxies. It is possible that our selection of E and SO galax- 
ies in the intermediate redshift clusters is biased to select 
already aged galaxies (see also Franx & van Dokkum 1996). 
In this case we are likely to underestimate the mean evo- 
lution measured from the FP when comparing to complete 
samples of galaxies in nearby clusters. Larger and more com- 
plete samples in intermediate redshift clusters are needed in 
order to address this problem in detail. 



6 CONCLUSIONS 

We have established the Fundamental Plane (FP) for the 
two rich clusters Abell 665 and Abell 2218, both at a red- 
shift of 0.18. The photometric parameters were derived from 
ground-based observations obtained with the Nordic Opti- 
cal Telescope, La Palma, and from HST observations. The 
photometry was calibrated to Gunn r in the rest frame of 
the clusters. Central velocity dispersions were measured for 
six galaxies in A665 and ten galaxies in A2218. The FPs for 
the two clusters were compared to the FP for nearby clus- 
ters derived by JFK96, and to the FP for the Coma cluster. 
The scatter around the FP for A665 and A2218 is similar to 
the scatter found for nearby clusters. 

We have used the data for A665 and A2218 together 
with other recently available data for intermediate redshift 
clusters (van Dokkum & Franx 1996; Kelson et al. 1997) to 
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test if the slope of the FP changes with redshift. The data 
for the clusters CL0024+16 (z=0.39), CL1358+62 (z=0.33) 
and MS2053-04 (z=0.58) were calibrated to Gunn r in the 
rest frame of the clusters. We find that the coefficient for 
the logCT term in the FP is significantly smaller for the in- 
termediate redshift clusters than for the Coma cluster. The 
difference cannot be explained by the difference in the lu- 
minosity distribution of the Coma cluster sample and the 
samples in the intermediate redshift clusters. However, other 
types of selection effects, like a difference in the distribu- 
tion of the combination log re -|- 0.82 log <7>o, have not 
been ruled out. The smaller value of the coefficient implies 



for the intermediate redshift clus- 
for nearby clusters 



-0.02 



that M/L oc M"-'"'ro 
ters, compared to M/L oc M'^'^'^r, 
(JFK96). This may indicate that the low mass galaxies in 
the clusters evolve faster than the more massive galaxies. We 
cannot distinguish if this is purely a difference in the lumi- 
nosity evolution or if also dynamical evolution takes place. 
Alternatively, the velocity dispersion and not the mass could 
be the driving parameter. Direct fits to the M/L ratio as a 
function of the velocity dispersion give M/L oc ct"'^® for the 
Coma cluster and M/L oc a^ *^ for the intermediate redshift 
clusters. 

The zero point offsets in the FP for the intermediate 
redshift clusters were used to investigate the average evolu- 
tion of the M/L ratios of the galaxies. The M/L ratios of 
the galaxies in A665 and A2218 are on average (11 ± 7)% 
smaller than the M/L ratios of the galaxies in the Coma 
cluster. From all the clusters we find that the M/L ratio 
changes with the redshift as A log M/Lr = (-0.26±0.06)A2 
(for go = 0.5), equivalent to a change in the absolute mag- 
nitudes of AMrT = (-0.65 ± 0.15)A2. The slow evolution 
of the M/L ratios as a function of redshifts that we find is 
in general agreement with previous results by van Dokkum 
& Franx (1996), Kelson et al. (1997), Bender et al. (1998) 
and Pahre et al. (1999). The result can be used to con- 
strain the formation redshift of the galaxies. The interpre- 
tation depends on the crucial assumption that the E and SO 
galaxies observed in the intermediate redshift clusters evolve 
into galaxies similar to present day E and SO galaxies (cf., 
Franx & van Dokkum 1996), and that the samples in the 
intermediate redshift clusters form a representative sample 
of all the galaxies that will end as present day E and SO 
galaxies. For go = 0.5, models with a single burst popula- 
tion, a Salpeter IMF and Zform > 5 are consistent within 2a 
with the observed evolution of the M/L ratios. Models with 
2form < 2.4 deviate more than 3cr from a least squares fit 
to the data. For go = 0.15, we find Zform ~ 2.5, while mod- 
els with Zform < 1-4 deviate more than 3(t from the fit to 
the data. The model constraints presented here are stronger 
than in previous studies due to the larger number of clusters 
included in the analysis. 
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Further, the low frequency variations over the fields were 
corrected for by using night time exposures of blank fields. 
Dust grains gave rise to flat field variations from night to 
night, and sometimes during the night. Night time expo- 
sures of blank fields were also used to correct for this effect. 
The pixel-to-pixel accuracy of the flat fields is 0.2-0.4% in 
areas not affected by the dust grains, and 0.8% in areas af- 
fected by the dust grains. The low frequency variations over 
the field and any residuals from the dust grains are below 
0.5% after the flat fleld correction. 

The I-band images of the two galaxy clusters were com- 
bined and cleaned for signal from cosmic-ray-events (CR) as 
follows. Geometrical transformations (shifts in x and y and 
a rotation) were used to transform all images of a field to 
a common reference image. The median sky level was sub- 
tracted from each image, and the relative intensity was de- 
termined from aperture photometry of objects brighter than 
~ 20™5 in the I-band. A median image was constructed, each 
image scaled according to the relative intensity. Each im- 
age was then compared to the median image and deviations 
larger than seven times the local noise were flagged as CRs. 
Mask images of the flagged pixels as well as other bad pix- 
els due to the CCD were constructed. Finally, a mean image 
was derived, omitting the flagged pixels, and each individual 
image scaled according to the relative intensity of the image. 
To maintain the statistically properties of the mean image, 
the mean of the median sky levels was added to the final 
image. The masks were inspected visually to make sure no 
signal from objects got masked. Further, the cleaned mean 
image wEis compared with a mean image with no masking 
by means of aperture photometry of the objects in the field. 

The magnitude zero points for the combined images 
were established by comparing aperture photometry from 
these images to aperture photometry derived from the in- 
dividual images. Photometry was derived for 40-60 objects 
with I-band magnitudes brighter than ~ 20'"5. The consis- 
tency of the zero points is better than 0^015. 



APPENDIX A: THE PHOTOMETRY 

Al Basic reductions of ground-based images 

The bias level and pattern turned out to depend on the CCD 
temperature, which in turn varied with ±1.0K. Therefore, 
the bias was subtracted row-by-row based on the overscan 
region of the images. After this processing a few of the im- 
ages showed row-to-row variations in the background of 1- 
1.5%. The level of these variations were determined directly 
from the images, and subtracted out with appropriate filter- 
ing to ensure the levels were not affected by signal from the 
objects. No correction for the dark current was performed. 
The dark current is less than 2e~ per hour. 

The offset between the set exposure time and the ac- 
tual exposure time for the iris shutter of the camera was 
determined from dome flats. The average shutter correction 
is A = 0.0826 ± 0.0008 sec. 

The linearity of the CCD was tested from dome flats. No 
signs of non-linearity can be detected within the detection 
accuracy of 1.5%, for levels up to the software saturation 
limit. 

Flat fields were constructed from twilight fiat fields. 



A2 Basic reduction of HST images 

The HST/WFPC2 images were corrected for bias, dark- 
current and flat fleld corrected using the pipeline reduction 
supplied by Space Telescope Science Institute. The images 
were co-added and cleaned for cosmic-ray-events using the 
iterative method implemented in the Space Telescope Sci- 
ence Data Analysis Software (STSDAS) task "crrej". Where 
necessary, the individual images were first registered by ap- 
plying small integer pixel shifts. 

A3 Calibration to standard passbands 

A3.1 Ground-based photometry 

The NOT photometry was standard calibrated based on ob- 
servations of stars in M67, M92, and the field SAlOl taken 
during the nights of the galaxy observations. All nights dur- 
ing which galaxy observations were obtained were photomet- 
ric. Aperture photometry was derived for 188 observations 
of a total of 35 stars. The observations have large variations 
in the seeing, and some of the exposures were taken with 
the telescope slightly defocused due to the brightness of the 
stars. The aperture radii for the photometry were typically 
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4'.'9. Based on growth curves of bright stars in the fields the 
aperture magnitudes were corrected for the finite sizes of the 
apertures. The growth curves typically extended to 9'.' 5, and 
the typical size of the correction is — CK"04. This correction 
ensures consistency between the defocused and focused ob- 
servations, and between observations with large differences 
in the seeing. All observations in the I-band were obtained 
during night f and night 5. All observations in the V-band 
were obtained during night 2. The typical seeing for the stan- 
dard star observations was FWHM=0'.'6-0'.'7, l'.'5-2'.'5 and 
0'.'7 — I'.'O for the nights 1, 2 and 5, respectively. 

The 9-15 stars in M67 were used for determination of 
the atmospheric extinction. We find A:i(nightl) — 0.063 ± 
0.001, ^/(nights) = 0.124±0.001, and fcv(night2) = 0.165± 
0.004. Further, observations from night 5 were offset to con- 
sistency with night 1 by subtraction of O^'Oll. 

We have used standard magnitudes from Landolt (1992) 
(SAlOl) and Davis (1994) (M92). The magnitudes from 
Davis are consistent with Christian et al. (1985) and Stet- 
son & Harris (1988) for the stars in common. The standard 
magnitudes from M67 are from Montgomery et al. (1993), 
Joner & Taylor (1990) and j0rgensen (1994), in this prior- 
ity. The magnitudes from the three sources are consistent 
for the stars in common. We derived the following standard 
transformations 



/in 



est 



0.021 



V = Vinst + 0.0596(y - /)i„st + est rms = 0.014 



(Al) 



The I magnitudes were calibrated to Cousins /c. The trans- 
formation for Ic has no significant color term; if we include 
a {V — I) term in the transformation we find a coefficient for 
this term of 0.001 ± 0.007. The photometry for the galax- 
ies was corrected for the galactic extinction. For A665, the 
galactic extinction is Ab ~ 0.15 equivalent to Ai = 0.07 and 
Av — 0.11. A2218 has a galactic extinction in the B-band 
of Ab = 0.10, equivalent to Ai = 0.05 and Ay = 0.08. The 
galactic extinctions were derived from the reddening maps 
published by Schlegel, Finkbeiner & Davis (1998). 



A3. 2 Photometry from HST observations 

The photometry from the HST images was calibrated using 
the transformations given by Holtzman et al. (1995b). For 
F814W we use the transformation based on observations, 
while for F606W and F702W we use the synthetic transfor- 
mations. We define instrumental magnitudes as 

minst = -2.51og(DN/fcxp)+ZP+2.51ogGRi+Am+0.05(A2) 

"DN" is the number of counts, "ZP" are the zero points 
given by Holtzman et al. ( 1995b) for the relevant transfor- 
mation (see equations |A4 and [A.5|). We use GR2 = 2.003, 



GR3 = 2.006 and GR4 = 1.955. The term 0.05 is added as 
a correction for the difference between "short" and "long" 
exposures (Hill et al. 1998). The aperture corrections Am 
were estimated from data for the encircled energy given 
by Hohzman et al. (1995a). We use Amgw = 0.0983, 
Am702 = 0.1087, and Ameoe ~ 0.1138. The instrumental 
magnitudes in the three filters are then given by 

771814 = -2.5 log(DN/fexp)-F 20.987 2.5 log GRi 
771702 = -2.5 log(DN/texp) -I- 21.670 2.5 log GR« (A3) 
777606 = -2.5 log(DNAexp) + 22.257 -1- 2.5 log GR« 



The photometry for galaxies in A665 was calibrated using 
the transformations for 7c and V and colors based on HST 
data. We use the following transformations 



(A4) 



Ic = 777814 - 0.062(V^ - Je) + 0.025(V - 
V = 777606 + 0.254(V - Ic) + 0.012(1/ - Jcf 

The Ic transformation is for the in-flight system based on 
stars bluer than the typical colors of the observed galaxies, 
{V — Ic) = 1.5. However, for {V — Ic) = 1-5 the transforma- 
tion agrees within 0™005 with the synthetic transformation. 

A2218 was observed in the F702W filter, only. We have 
calibrated this photometry to Ic using the transformation 
for Rc given by Holtzman et al. (1995b) and the relation 
{V — Rc) ~ 0.52(V — Ic), which we derived from standard 
star magnitudes for stars with (V — Ic) in the interval 1-2 
(Landolt 1992). The coefficient is in agreement with colors 
of E and SO galaxies at 2; = 0.2 (e.g., Fukugita, Shimasaku 
& Ichikawa 1995; Frei & Gunn 1994). This gives the trans- 
formation 



771702 - 0.227(1/ - Ic) - 0.021(1/ - Ic)'^ 



(A5) 



We used colors from the ground-based photometry in order 
to calibrate the HST observations of A2218. Finally, the 
photometry was corrected for the galactic extinction. 

A4 Photometric parameters 



Tables Al and A2 give the effective parameters and colors 
of all galaxies in the two clusters brighter than 777t = 19™ in 
Ic and within the fields covered by the NOT observations. 
Five of the galaxies in A665 were also observed by Oegerle 
et al. (1991). The cross-references are as follows, with our 
galaxy number given first. 1032 = 218, 1150 = 201, 1168 = 
224, 2105 = 225 and 2122 = 227. 



A5 Comparison of ground-based and HST 
photometry 

The ground-based photometry and the HST photometry 
were compared by means of aperture magnitudes. In or- 
der to limit the effect of the different spatial resolution and 
the choice of the aperture size we convolved the HST im- 
ages to approximately the same resolution as the ground- 
based images. Aperture magnitudes for the galaxies were 
then derived within apertures of with radii of 3 arcsec. The 
magnitudes wer e sta ndard calibrated as described in Sec- 
tions A3.1 and A3. 2, Fig. Al show the comparison of the 



aperture magnitudes. We find the mean differences to be 
/c(NOT) - 7c(HST) = 0.066 ± 0.007 and 0.017 ± 0.010 for 
A665 and A2218, respectively. Galaxies with aperture mag- 
nitude fainter than 19™ were excluded from this comparison 
in order to limit the errors introduced by uncertainty in the 
sky subtraction. 

The effective parameters derived from the NOT images 
and the HST images were compared, see Fig. and Table 



A3, Only the galaxies with total magnitudes brighter than 
19™ in the 1-band were included in this comparison, a total 
of 56 galaxies. For the fainter galaxies the individual param- 
eters (logTc, <M>c and 771t) become highly uncertain due to 
the small size of the galaxies compared to the FWHM of the 
PSF for the ground-based data, and due to rather low signal- 
to-noise for these faint galaxies in the HST observations of 
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TABLE Al 

Photometric parametkrs, Abell 665 



R.A. 
(J20Q0) 



Dec. 
(J2000) 



Field 



log re 
farcsecl 



— NOT 

<M>e 



(V - h) 



log re 
farcsecl 



<M>e 

(F606W) 



(V) (F606W) 



- HST ■ 
log re 



<M>e 

(F814W) 



<M>e 



(F814W) 



(V - /c) 
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31 


00 
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65 


49 


00 


2 


1 


-0 


248 


19 


30 


1 
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49 


-0 


172 


20 


71 


21 


01 


2.3 


-0 
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19 


63 


19 


59 


1 


9 


1.46 


1032 8 


31 


10 


37 


65 


49 


17 


3 


1 


-0 


010 


19 


15 


3 
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1 
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06 
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49 
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1 





091 


19 


92 


1 


2 


1 
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NOTE. — Galaxy identifications - this paper. Field - the NOT field that contains the galaxy, see Table 3 for seeing values. <fj,>^{Ic], <i^>q{V) and {V — Ic) for both the NOT and the 
HST data are standard calibrated to the Johnson-Kron-Cousins photometric system and corrected for galactic extinction. <M>e (F814W) and <fS>^ (F606W) are the instrumental magnitudes 
in F814W and F606W, respectively, see Eq. (A3). The formal fitting uncertainties are typically as follows, log ±0.006 for NOT data, ±0.002 for HST data. <(x>^: ±0.02 for NOT 
data, ±0.007 for HST data. The typical (random) uncertainty on the {V — Ic) is 0.01. The goodness-of-fit is given as x^/^pix, where A^pix is the number of pixels within the fitting radius. 

— ^(*^« ~ ^modex)^ / f^f with the sum over all the pixels within the fitting radius. Ui is the signal in pixel «, nj^modei is the model value in the pixel, and erf the noise (photon- and 
read-out-noise) in the pixel. All the galaxies have been fitted with r^^^ models. 
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TABLE A2 
Photometric parameters, Abell 2218 



NOT HST 

Galaxy R.A. Dec. Field log rc <ii>^ X^/Np-,^ {V - Ic) logr, <ij->^ <At>e X^/JVpix 

(J2000) (J2000) [arcsec] (/e) [arcsec] (F702W) (/„) 
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Note. — Galaxy identifications - Le Borgne et al. (1992) (L numbers), Butcher et al. (1983) (BOW numbers), this paper (J number). 
Field - the NOT field that contains the galaxy, see Table 3 for seeing values. <;j>^(/<;) and (V — Ic) for both the NOT and the HST 
data are standard calibrated to the Johnson-Kron-Cousins photometric system and corrected for galactic extinction. <ii>^ (F702W) is 
the instrumental magnitude in F702W, see Eq. (A3). The formal fitting uncertainties are typically as follows, log rc: ±0.006 for NOT 
data, ±0.002 for HST data. <IJ.>^: ±0.02 for NOT data, ±0.004 for HST data. The typical (random) uncertainty on the (V - /„) is 
0.01. The goodness-of-fit is given as X^/^pix, where JVpix is the number of pixels within the fitting radius, = J^("i ^ , model ) ^ /f,^ i 
with the sura over all the pixels within the fitting radius, nt is the signal in pixel i, n;, model is the raodel value in the pixel, and trf the 
noise (photon- and read-out-noise) in the pixel. The galaxies L114, L395, BOW83 and BOW97 have been fitted with exponential disk 
models, all other galaxies have been fitted with r^'^ models. 
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Figure Al. Comparison of aperture magnitudes. An aperture 
with radius 3 arcscc was used. 



Table A3. Comparison of effective parameters 



Parameter 


mean difference 


rms 


log rc 


-0.033 ± O.OfS 


O.lf 




-0.fl9 ± 0.053 


0.4f 


mx 


0.044 ±0.02f 


0.f6 


ppa 


0.006 ± 0.003 


0.026 



Notes. " the combination logrc — 0.328</i>(,. The differences are 
calculated as "NOT"-"HST". Only galaxies with total magni- 
tudes brighter than f9™ were included in the comparisons. The 
two clusters are combined in this table, and the fIST photometry 
has been offset to consistency with the ground-based photometry, 
see text. 



A665. The combination "FP"=logrc - 0.328</i>c, which 
enters the FP, was also compared. This combination is well 
determined even for the fainter galaxies, but since none of 
the galaxies for which we have spectroscopy are fainter than 
19™ we restrict the comparisons to galaxies brighter than 
this limit. The rms scatter for the comparison of the combi- 
nation "FP" is similar to typical values for comparisons of 
data for nearby clusters (e.g., j0rgensen et al. 1995a). If the 
uncertainty is equally distributed on two data sets, the typi- 
cal uncertainty on the combination "FP" is 0.018. The small 
mean difference in the combination "FP" between the NOT 
data and the HST data is equivalent to a difference in the 
FP zero point of only 1.4%. The comparisons of the param- 
eters logrc, <M>c ^ii'i "^T have slightly larger scatter than 
usually found for comparisons between effective parameters 
for nearby galaxies derived by different authors. 

The errors on logrc and <At>c s^re highly correlated 
(e.g., j0rgensen et al. 1995a). We use this fact to provide 
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Figure A2. Comparison of effective parameters in the I-band. 
Boxes — A665; triangles - A2218. Solid symbols — galaxies for 
which we have spectroscopy. The differences are calculated as 
"NOT"-"HST". The formal errors derived from the fitting are of 
the size of the points. "FP" is the combination log rc — 0.328<^>c 
which enters the FP. Only galaxies with total magnitudes brighter 
than 19™ are shown. The HST photometry shown on the figure 
has not been offset to consistency with the ground-based photom- 
etry 



another chec k of the magnitude zero points of the HST pho- 
tometry. Fig. A3 shows the difference in log rc versus the dif- 
ference in <fj,>^. The differences were calculated as "NOT"- 
"HST". The dashed lines on the panels show least squares 



fits to the data. If the zero point for the HST photometry 
was in perfect agreement with the ground-based zero point, 
then the dashed lines would go through (0,0). The offset in 
A</i>^ is equal to the difference in the zero points. We find 
7c(NOT) - /c(HST) = 0.063 ± 0.011 and 0.015 ± 0.008 for 
A665 and A2218, respectively. This is in agreement with the 
determinations of the zero point differences based on aper- 
ture magnitudes. 

Offsets of a few hundredths of a magnitude are com- 
mon between HST photometry and ground-based photome- 
try (e.g., Ellis et al. 1997; HiU et al. 1998). The offset we 
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Figure A3. The difference in logrc versus tlie difference in <fi>^. 
The differences are calculated as "NOT"-"HST". Only galaxies 
with total magnitudes brighter than 19™ are shown. Solid symbols 
— galaxies for which we have spectroscopy. The dashed lines show 
least squares fits to the data. The offset in A<fj,>^ can be used to 
derive the offset between the HST photometry and the ground- 
based photometry, see text. 



find for tfie A665 pliotometry is sliglitly larger tlian ex- 
pected. Ifowever, in order to ensure tlie consistency of the 
data for A665 and A2218 we have chosen to caUbrate the 
HST photometry to the ground-based /c- We have applied 
the following offsets: h = /c(HST) + 0.065 for A665 and 
Ic = /c(HST) + 0.016 for A2218. 

We have no reliable way of checking the zero point 
for the V magnitudes based on the data, because of 

the poor seeing of our ground-based V-images. We have 
therefore not applied any zero point correction to the colors 
{V - Ic) derived from the HST data. 



APPENDIX B: SPECTROSCOPY 

The spectroscopic observations were made through multi- 
slit masks. One field was observed in each cluster. The final 

and 



spectra are shown in Figs. Bl 



B2 



The reductions of the 
A2218 spectra are described in the following. The methods 
for the reductions of the A665 spectra were very similar. We 
refer to Franx (1993ab) for additional information regarding 
the A665 spectra. 



Bl A2218 

The spectra were corrected for the bias and the dark by 
standard methods. The bias subtraction was based on the 
level in the overscan region as well as individual bias frames. 
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Figure Bl. Spectra of the observed galaxies in A665. The spectra 
are not flux calibrated. 



The flat fields show strong variations in the fringing 
depending on the zenith distance of the telescope due to the 
flexure of the instrument Therefore, flat flelds obtained in 
day time did not give sufficiently good flat field correction. 
Instead we use flat flelds obtained during the night. These 
flat fields were obtained with the internal quartz lamp in the 
instrument. Two flat fields bracket each science exposure. 
The flat fields obtained in day time were combined in the sets 
they were taken cleaned for cosmic-ray-events (CR). Then 
all the flat fields were normalized by fitting a high order 
spline in the direction of the dispersion. In order to clean the 
night time flat flelds for CRs they were each compared to 
the the day time flat field, which best resembled night time 
flat field. The final flat flelds were derived as the average of 
the two fiat fields, which bracket the science exposure. The 
typical pixel-to-pixel noise in the final flat fields was 0.6%. 
Areas of the flat flelds that had a pixel-to-pixel noise larger 
than 1.5% were set to one, and no flat field correction was 
performed. This also ensures that no artificial variations are 
introduced at the edges of the slit-lets. 

The individual exposures of A2218 were cleaned for CRs 
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as follows. First, the median image was derived. Then differ- 
ence between each individual image and the median image 
was calculated. This image naturally has strong residuals 
from the sky lines and the galaxies. The residuals from the 
galaxies wore fitted with a smooth function in the direction 
of the dispersion, and the fit was subtracted. The residuals 
from the sky lines were fitted with a constant across each 
slit-let, and the fit was subtracted. The resulting difference 
image contains noise and signal from CRs, plus some resid- 
uals from the strongest sky lines. Pixels for which the devia- 
tions were below 7 times the local noise or which had resid- 
uals from the sky lines were set to zero. This resulted in a 
mask image for each of the science exposures. The mask im- 
age was subtracted from its corresponding science exposure. 
The net-effect is that pixels affected by the CRs get substi- 
tuted with the median value, variations due to the shifts of 
the galaxy spectra and the sky lines taken into account. 

The spectra were corrected for the slit function based 
on sky flat fields. The derived slit function was shifted up to 
0.4 pixel in the spatial direction in order to take into account 
the shifts in the slit-lets due to the flexure. 

The wavelength calibration was done for each slit-let in- 
dividually. The main wavelength calibration was established 
from high signal-to-noise He-Ne-Ar lamp spectra taken in 
day time. The dispersion solution was then shifted slightly 
based on He-Ne-Ar lamp spectra taken immediately before 
and after each science exposure. The typical rms of the wave- 
length calibration is O.ISA. 

The spectra were corrected for the S-distortion by trac- 
ing the position of the spectrum of the BCG. Because the 
signal in the individual exposures is too low to establish the 
S-distortion for the fainter galaxies the same S-distortion 
was used to correct all the galaxies. After the correction was 
applied, we verifled that the same correction could be use 
for all the galaxies by deriving the differences between the 
mean center of the the spectra in the wavelength intervals 
5295-5445A and 6040-6130A. The BCG has a difference of 
0.02 pixels, while the rest of the galaxies show differences 
of less than 0.2 pixels. The aperture used for the velocity 
dispersion measurements is 5 pixels. We therefore conclude 
that the adopted correction for the S-distortion does not 
affect the velocity dispersion measurements. 

The wavelength calibration and the correction for the 
S-distortion were applied simultaneously resulting in indi- 
vidual science exposures that are all on the same wavelength 
scale and with the galaxies centered on the same pixels. The 
spectra were sampled on a logarithmic wavelength scale. The 
individual images were then averaged with scaling according 
to their exposure time. The individual expsure times were 
11 X Ih and one exposure of 48 minutes. The average im- 
age was at this stage cut into sub-images, each sub-image 
contains the spectrum from one slit-let. The sub-images are 
then treated as individual long-slit spectra. 

The sky background was subtracted from each sub- 
image based on the outer 3'.'5-10'.'35. When possible the sky 
spectra were extracted from both sides of the galaxy spectra. 
In some cases residuals were left from the stronger sky lines. 
The sky subtraction of these lines was improved by fitting a 
linear function to the residuals in the spatial direction, and 
then subtracting the fit. The residuals from the sky line at 
5577A were interpolated across. 

The spectra were calibrated to a relative flux scale based 



on exposures of the flux standard star HD192281. The flux 
calibrated spectra are in this paper used for display purposes 
only. The determination of the velocity dispersions does not 
depend on the flux calibration. 
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Figure B2. Spectra of the observed galaxies in A2218. The spectra are calibrated to a relative flux scale. The panels are labeled with 
galaxy numbers from Le Borgne et al. (1992) or Butcher et al. (1983). BOW97 is a background galaxy at z=0.291. 
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